x1=[0.7	0.8	1	1.1	1.2	2.81	3.8	2.99	2.1	1.5	0.9	0.57	0.3	0.42	0.5	0.41	0.4	0.69	0.9	0.86	0.8	0.74	0.7	0.79	0.9	1.14	1.2	0.76	0.4	0.38	0.4	0.57	0.7	0.74	0.7	0.42	0.3	0.72	1	0.84	0.7	0.7	0.7	0.8	0.8	0.58	0.4	0.4	0.4]';
x2=[0.7085	0.814	0.9551	1.1879	1.2229	1.5118	3.8054	3.2577	2.2999	1.5567	1.1679	0.4842	0.5304	0.0742	0.6283	0.3147	0.4161	0.3296	0.8676	0.7828	0.78	0.6881	0.6637	0.6194	0.8007	0.8442	1.2096	1.0112	0.4437	0.3041	0.3818	0.3542	0.6529	0.6331	0.7166	0.581	0.2287	0.3121	0.9403	0.9036	0.6681	0.6252	0.6608	0.6348	0.8221	0.6721	0.4262	0.3238	0.4052]';
x3=[0.7166	0.8135	1.0159	1.1134	1.3206	1.7276	1.7085	3.3141	2.6653	1.9194	1.2771	0.9814	0.4162	0.4181	0.0673	0.5126	0.2613	0.3378	0.2753	0.7232	0.6606	0.6511	0.5763	0.5585	0.5288	0.6871	0.7264	1.0216	0.8565	0.3756	0.257	0.3249	0.3023	0.5584	0.5397	0.6073	0.4969	0.1955	0.2688	0.8	0.7725	0.5743	0.5381	0.5724	0.5485	0.7069	0.5784	0.3674	0.279]';
x4=[0.4636	0.7933	1.0918	1.1682	1.1608	1.9524	2.13	0.6179	0.7207	2.1762	1.7334	1.3041	0.8916	0.6517	0.3011	0.2798	0.0501	0.3523	0.1937	0.2374	0.1981	0.5038	0.4647	0.4666	0.4252	0.4155	0.3966	0.4892	0.5091	0.7324	0.6114	0.2731	0.1887	0.2375	0.2208	0.4001	0.3902	0.4462	0.3697	0.1425	0.1957	0.5947	0.5758	0.4334	0.4048	0.422	0.4039	0.5256	0.4304]';
num=length(x1);
res=0
for i=1:20
   res=res+(abs(x4(i,1)-x1(i,1))/x1(i,1));
end
res/40
plot(1:num,x1,'r+-',1:num,x2,'b+:',1:num,x3,'bo:',1:num,x4,'bx:');